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ABSTRACT 


Three monthly mean oimulationo of the global atmoo- 
phere were computed for February 1976 with the GISS model 
from obaerved Initial conditions on the first day of the 
month. In a replication experiment, two of these com- 
putations generated slightly different monthly mean 
states, apparently duo to the schedule of interruptions 
on the computer. The root-mcan-square errors of replica- 
tion over tli^ Northern Hemisphere were found t*^ be about 
2 mb, 20 m, and 1 K for sea-level pressure, 500 mb height, 
and 850 mb temperature, respectively. The monthly mean 
500 mb forecast results for February 1976 over the Northern 
Hemisphere are consistent with those from earlier GISS 
model experiments, and again indicate some predictive 
skill at that level. Use of the observed monthly mean 
sea-surface temperature (SST) field for February 1976 in 
place of the climatological SST field for February re- 
sulted in slightly improved simulations over the globe 
cuid Northern Hemisphere, but not over smaller subregions. 


1 



I 


1. Introduction 




Mils study is a pdiitinuation of an effort to deter- 
mine if a global general circulation model (GCM) initial- 

r, '' 

ized v/itb data for the first day of a month can simulate 
realistically the mean state of the atmosphere for that 
month. Eaiiler experiments in monthly mean prediction 

ii f’’ ■ /) 

v>, jl ’ 

(Spar et al., 1976; Spar, 1977) with the GCM kno^m as 
the GISS^ model (Somerville et al., 1974) have ^indicated 
that, despite the well-lmown decay of predictability, 
time-averaging of a forecast history over a month may 
rcsu.Vb in some modest improvement over climatology. How- 

ever;' that preliiiiinary assessment was based on only tliree 

i’i 

cases: January of 1975, 197 'V and 1975. In this note 

we present the results of still another set of monthly 
mean simulations with the model, for February 1976. 

As in the earlier computations, the Pebruary 1976- 

7. S' 

global data set used both for initialization (00 GMT, 

1 Pebruary) -and verification was provided by the National 
Meteorological Center (NI4C) from its operational 
archive. A climatological mean Pebruary atmospheicic 

global data set employed in part of the forecast evalua- 

" 7. i .-“ 7 

tion was obtained from thi^N^ Center for Atmos- 

pheric Research (NCAR), which also furnished the ’ 


^The GISS model was developed at the Goddard Institute , 
for Space Studies. 
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climatological mean February sea-surface temperature (SST) 
field (V/ashington and lEhiel, 1970) used in two of the 
three separate February 1976 simulations,, ‘^In the third 
computation, a monthly average SS5D field for February 
1976, based on satellite scanning radiometer data pro- 
vided by the National Environmental Satellite Service 

fi 

(NESS), replaced the climatological SSI’ field, 

QIhe first two simulations, identified below as FI 

and F2, were computed as part of a replication expert- 

- '''' ?;\\ 

"1 ■ 

ri ment. Run on the same IJ3M 360/95 computer with identi- 

■i 

cal programs, boundary conditions, and initial condi- 
“tions, FI and F2 presiMably differ only because of dif- 
ferences in the schedule of interruptions and restsirts 
necessitated by other -demands on the computer during 
the month-long forecast.' A comparison of forecasts FI 
and F2 thus provides . a measure of the reproducibility 
of results, IDhe third simulation, design^, ted F3, gives 
some indication,;/ when compared with FI ,ef F2, of the in- 
fluence of ocedt’A surface temperature variations on the 
predicted monthly mean atmospheric fields. 

As discussed by Dickson (1976), .February 1976 v/as 
characterised by unusual, ly warm weather over most of the 
United States, v/here new February high temperature records 
- were set at numerous locations, notably > in the east and 
southwest, Ihis high temxaerature anomaly v/as associated 


Qli' POOR QUAUTXI 


with the development during the month of a persistent, 
fast, flat zonal tropospheric flow pattern in middle and 
high latitudes over the western half of the Northern 

n 

Hemisphere, accompanied by a northv/ard shift of the ' 
strongest westerlies. Only over eastern Europe and a 
westem Asia did the hemispheric mean tropospheric cir- 
culation for the month exhibit the large amplitude short 
v/ave patteir normally associated with strong meridional 

thermal adVi '.:tion. Over North America', on^the other 

' '' 

hand, a straight westerly current of mild{ maritime air 
_ from the Pacific, per sis tod f^^ most of the month, with 
little southv/ard transport of cold continental air from 
Canada, As a result, temperatures remained abnormally _ 
high during the month over much of North America. [Hiis 
monthly mean circulation pattern and the associated tem- 
perature anomaly were not In place at the "beginning of 
the month, but developed after the first- week. IDhus, 
February 1976 is a challenging test of the ability of 
the model to simulate - a major circulation and temperat-ure 

anomaly from initial conditions in v/hich the monthly mean 

\\ . ' ' ^ ■ 

' pattrern was not yet apparent. 

2. Replication Experiment 

As in ovcc earlier forecast studies, the results/ of 
the February 1976 simulation experiments are presented; ■ 
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■belqtf in terms of onergetics, mean wind profiles, fore- | 
cast error statistics, and monthly meaH'- synoptic maps, / 
with emphasis on the Northoi’iinifemispliers,. 

In the replication experiment, the two monthly moan 
simulations for Fehruar*:^ 1976, designated PI and P2, were 

carried out with the same program, and with identioal 

' f 

initial and houndgiry conditions, on the same coisipiiter. 

One might therefore expect identical results i^rom these 
tv/o calculations. Hdv;ever, as it v/as not possihle to 
carry out such a long simulation run in one continuous 
operation at GISS, because of other demands on the com- 
puter, the month-long forecasts were computed piecemeal, 


with several interruptions and restarts. (Uhe interrup- 
tion schedule v/as .indeed different for PI and P2, and 
this may account for the differences between the two re- 
sults. Because of the cycling schedule in the model 


program, and the fact that not 7all calculations 

■jL 

identical, in each time step (Somerville et af. , 1974) > 
it is apparently possible for small differences to de- 
velop betv/een tv/o long, interrupted forecast runs. She 
purpose of this; experiment was to determine the magni- 
tude of these differences. 

In Table 1 are shown the forecast (Pl and P2) and 
observed (O) global, and northern hemisphere mean ener- 
getics for Pebruary 1976. and Kjyj represent the mean 


ga~blo 1 . Zonal availalsle potential energy zonal 

kinetic energy and eddy available 

potential energy (Pj,) and eddy kinetic energy 
(Kjj) of standing waves only, for Pobruai'y 1976 
over ,the Northern Hemisphere and the„_globe, 
integrated up to the 120 mb level 'for fore- 
casts PI and P2 and the observed (O) mean 
monthly atmosphere. Units: lO-' J m . 


Korthern Hemisphere ■ G-lobe 



' PI 

P2 

0 

PI 

P2 

0 

% 

58.9 

58.0 

49.8 

42^5 

41.8 

55.9' 

% 

8.56 

, 8.42 

7.21' 

6.69 

0.61 

5.97 


2.55 

2.45 

_5^71 

1.69 

1.71 

2.29 


1.14 

1.51 

2.25 

0.91 

1.00 

1.54 
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zonal available potential energy and mean zonal kinetic 
energy, respectively, integrated up to aijproximately the 
120 ml) level (the base of the highest model layer), while 
Pj, and ICj, represent the corresponding eddy energies for 
the montlily mean staiiding waves only, all in units of 
10^ J ra“ . (Por further details on the model energetics, 
see Somerville et al,, 1974; 'i'enenbaura, 1976; Spar et al. 
1976; Spar, 1977). 

,'V 

Prom fable 1 it is clear that the forecast mean,, 
energies are not identical in the two simulations. Hov/- 

ii 

ever, the differences between PI and P2 are small rela- 
tive to the ’'errors" of either forecast as measured 

l"! 

against the "observed" atmosphere. Hie relative dis- 
crepancy between the corresponding energies of PI and 
P2 is only about 10^ of the errvcr,, with the maximum dif- 

‘ / t \ i'' 

ference (17/“ of the error) in the eddy kinetic energy 

■ '' ^ , ■ / 

over the Northern Hemisphere. ' 

I 

I As in the previous monthly mean forecasts (Spas et 
al,', 1976; Spar, 1977), the G-ISS model overpredicts the 
zonal mean kinetic and potential- energies and under- 
predicts the eddy energy, fhese results a^e not unex- 
pected, as the model is ]movm. to generate Arctic i-zinter 
te-iperatures that, are too low due to its underestimate 
of- the eddy heat transport to high latitudes. The re- 
sulting overestimate of the meridional temperature 


7 


gradient produces excessive aoixE^ available potential 

.1 i‘ 

energy, as well as excessive westerly winds and aonal 
kinetic energy, IHio \mder estimate of eddy energies isi 
characteristic of all coarse-grid G0M*s. (V/ith j.’espect 
to the eddy energies, the ,rc»mlts for ITebruary 1976 are 
not strictly comparable to those previously published, 
as the latter included the 12-hoxa'ly transient eddy 
energies averaged over the month, while in bhe present 
study only the eddy energies of the monthly mean state, 
i.e. , of the standing eddies, are considered.) 

fhe predicted (3?! and P2) mean meridional profiles 
of the vertically- integrated aonal wind for the nine 
model layers up to thenlO mb level are showi in 
Tigure i, together with the observed profile based on 
the February 1976 MC global analyses. Wiile the -two 
simulations lie very close to each other, they do differ 
in some details, fhus, l^l indicates a mean v/ind maximum 
of 20,1 ins"*^ at latitude 5(|‘° H while 3?2 shows a maximum 
of 19.7 ino”^ at latitude , iEhe observed mean wind 

profile for tlie liorthern Hemisjphere shows that* both .com- 
putations seriously ovex'pf edict '^the maximum mean wester- 
lies in middle latitudes, wliich reached only 1^.1 ms"^ 
at latitude 58° N, whileamderprCdicting the westerlies 
in higher latitudes. In the Southeim Hemisphere both 
forecasts of the majcimum mean westerlies lie closer to ‘ 
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iig. 1. Mean taeridional profiles of vertically averaged r.eaR ror-al winds for February* 197£. 
FI and F2 represent the n:odel sirtulations, and 0 is the observed profile. 



obaervud value, liilth are ahifted one grid length 
(4°) too iar north. 

Both the Bl and P2 siimilationa generate a v/est- v/ind 
inaximiua at -the nominal 175 mh level oi the model at 
latitude 50° N in the mean meridional sonal, v/ind cross- 
section (not sho\m) . I'he location of the predicted suh- 
tropical jet stre^ is in good agreement v/ith the ob- 
served jet, which is also at. 50° N. However, the pro- 

•^*1 ~ 

dieted v/ind maxima are 36.5 tES for PI and 55.5;" Is 
for P2, while the observed mean maj:imum is 35»1 ms . 
tdie differences between the two simulations, as v/ell as 
their deviations from the observed v/ind profile, are 
illustrated in'Pigia’e 2, v/hich shows the tlireo vertica3,' 

profiles at latitude 50° H. Hers- it can be seen that 

-1 

forecast differences greater than 1 ms in the monthly; 
mean sonal v/inds may res,^t purely from computational 
discrepancies. The observed wind profile in Pigure 2 
further illustrates the fact that both forecast computa- 
tions overpredicted the mean westerlies at all levels, 
although the error at the jet stream level is q.uite 
small. _ , 

'JKie disci'epancies between forecasts PI onl P2 may 
also be expressed in terms of the root-mean-sq.uar’e (rms) 
differences and SI comparison scores (i’eweles and V^obus, 
1954) for various fields. (Poi' both statistics, perfect 


identity of fielda is ropreoentod "by values of sero. How- 
ever, ajx SI score of 20 io commonly rated as “virtually 
perfect" by MO forecasters, vdiile an Si score greater 
than 70 is generally considered “v/ortliloss,") In lable 2 
these statistics are shown for seven regions of the 
earth; (l) the globe; (2) the Northern Ilejuisphore ; 

(3) the “tropics" between latitudes 22® N and 22® S; 

(4) an East Pacific-Hnited States band (latitudes 30® N- 
54® N, longitudes 75^ V^-180°); (5) “North America" (lati- 
tudes 30® N-7<"'' longitudes 75° W-130® V/); (6) the 
"United StaieS" (latitudes 30° N-54® N, longitudes 75® Vi'- 
130® N); and (7) "Eumpe" (latitudes 54® N-S6® N, longi- 
tudes 10® VJ-40® E) . Ihe rms differences and SI compari- * 
son scores for El vs. 12 are shown in ffiable 2 for sea- 
level pressure and 500 mb height for all seven regions 
(with the omission of SI scores for region 4) ^ The rms 
differences of 850 mb temperatiires are shovm only for 
regions 2, 4, and 6, 

fable 2 provides a provisional estimate of the mini- 
mal error in a inontlily mean forecast that can be attri- 
buted to computational uncertainty alone, i . e , , the 
error res\;lting purely from the inability of the comput-er 
system to replicate perfectly a monthly me-an simulation 
from identical initial, and boundary conditions, (It 
does not, of course, reflect that additional computational 
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(Dable 2. 



Root-inean-square (rms) differences and SI \ 

comparison scores for Pebruary 1976 fore- 
casts PI vs. P2. 

jj /•, 



■'1 \ Re/ 2 :ion 

Sea-level 

Pressure 

500 mb 
IIei/3Cht 


850 mb 
Temperature 

rms (mb) 

81 

rms (m) 


rms (IC) 

1. ’ 

Globe 

1,5 

58 

16 

21 


2. 

Northern Hemisphere 

1.8 

43 

19 

23 

1.2 

5. 

I’ropics 

0.8 

37 

5 

33 



E. Pacific~U. S, 

2.1 


28 


1.2 

5., 

North America 

. 2.5 

51 

26 

20 


6. 

tmited States 

2.2 

48 

28 

16 

- 1.4 

7. 

Europe 

2.1 

43 

19 

25 

- " 


tmcertalnty associated with the existence oi alternative 
numerical approximations.) For sea-level pressure, 500 mb 
height, and 850 mb temperature, these minimal rms errors 
are approximately 2 mb, 20 m, and 1 K ior the Northern 
Hemisphere, while the minimal SI scores for sea-level 
pressure aiid 500 mb height are approximately 40 and 20, 
respectively, over the hemisphere. For the present, 
these may be regarded as the irreducible errors in the 
raontlily mean forecasts computed v/ith the GISS model. 

0 

5. Forecast Verification Results 

(Che forecast verification scores for FI and F2, 
based on the "observed" monthly mean state for February 
1976 as derived from 12 -hourly operational NHC analyses, 
arc listed in I’able 3. Also shown are the scores for a 
"forecast" of "Climatology" (M), based on a climatologi- 
cal mean February data set provided by NCAR. !The veri- 
fication regions are the .same as in (Dable 2, 

As might have been anticipated from Table 1, the 
errors of the FI and F2 simulations are similar, but not 
identical. Over the whole Northern Hemisphere, the dif- 
ferences between FI and F2 scores are negligible. Thus, I 
for example, the hemispheric rms errors in sea-level 

^ j . 

pressure sand 850 mb temperature are the same (8.8 mb and 
4o4 K) for_J5oth, and the hemispheric rms errors in 500 mb 
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Table ? 


, Hoot-mean-oquare (rms) erroro and SI oklll 
for February 1976 forccar>ta PI and P2. Aloo 
shovm are the ocorcs for forecaats < T 
Climatology (M). (See text for details.) 


A. Sca~lcvcl Preaoure 




n 


F? 


14 



Region 

rmM (mt>; 


rma (mb) 

Si 

rma (mol 

si 

1. 

Globe 

7.5 

66 

7.6 

67 

5.6 

71 

2. 

Northern Hemiophero 

8.8 

75 

0.0 

76 

6.3 

01 

3. 

Tropica 

5.7 

58 

5.8 

57 

2.7 

71 

4. 

E. Pacific-U. S. 

5.0 


4.9 


4.3 


5. 

North America 

7.3 

86 

7.0 

04 

5.4 

105 

6. 

United Statea 

5.0 

94 

5.0 

88 

4.5 

105 

7. 

Europe 

18.0 

113 

18.3 

117 

15.3 

117 


B. 

500 mb 

Hei/'ht 







PI 


F2 


M 



Region 

3*ins (nj 


rma (m) 

To*" 

r.”’.; (ill) 

S.L 

1. 

Globe 

66 

39 

60 

41 

79 

49 

2. 

Northern Hemisphere 

80 

42 

02 

44 

99 

53 

3. 

Tropica 

32 

59 

33 

58 

35 

60 

4. 

E. PacirJe-U. S. 

101 


110 


90 


5. 

North i\mcrica 

119 

37 

125 

38 

97 

51 

6. 

United States 

127 

20 

133 

29 

90 

48 

7. 

Europe 

109 

72 

191 

75 

234 

79 


C. 

850 mb 

Temperature (rma 

1 error: K) 



Region 

FI 


P2 


M 


2. 

Northern llcmiephcre 

4.4 


4.4 


4.6 


4. 

E. Pacific-U. S. 

5.3 


5.7 


4.6 


6. 

United States 

6.1 


6.4 


5.6 
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height (80 and 82 ai) differ only slightly. !Clie same is 
also true of the SI scores. Hov/ever, over smaller regions, 
such as the United States, the differences between the 

errors of the two simulations are -somtiv^at larger. 

/! ' ; ' 

!me rras differences and SI comparison scores of PI 
vs. P2 shovni in fable 2 may now be compare^wi-tU- the 
corresponding foi^ecast errors in fable 3. Over the 

' li " " - 

Northern Hemisphqre, for example, the rms difference of 
sea-level pressure between PI ^d P2 is 1,8 mb compared 
with an rms forecast error of 8,8 mb, while for 500 mb 
heights the corresponding values are 19 m and 81 ra.“ fhus, 
in both cases the rms replication error is approximately 
20J^ of the rms simulation error, Por 8501,ib temperature, 
the replication and simulation errors over ,’^he Northern 
Hemisphere are 1.2 and 4»4 K, respectively, and the rat'io 
is about 27!^. Although it is certainly not trivial, the 
replication problem may be regarded as minor compared 
with the prediction problem, at least in terms of rms 
errors, fhe problem of replication appears somewhat more 
serious vdth respect to the reproduction of gradients, 
particularly in the case of the sea-level pressure field, 
where the SI comparison score between PI and P2 over the 
Northern Hemisphere .Us 43 compared with SI simiilation 
scores of about 75. At 500 mb the hemispheric SI spore 

ii 

for replication is also more than hslf the magnitude of ' 
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the SI simulation score, hut in this case the value oi the 

'' 

replication score (25) is small enough to be regarcted_as 
negligible. " ^ ^ 

Sttie ’’Climatology" forecast (M) provides a standard 
against which the model simulations may be measured. As 
shov/n in lable 5, the model simulation for February 1976 
is not "consistently superior to climatology. In terras 

of rras error, climatology provides a better simulation's . 

il 

of the February 1976 sea-level pressure field over both 
the Northern Hemisphere and the globe than does the model. 
On^ the other hand, at 300 mb the model simulation over 
both the hemisphere^ and the globe , is superior to climC- 
tolo^ in terms of both rras errors and SI shill— scores, 
although this positive result is not found consistently 
over smaller regions. 

lo put the simulation statistics for February 1976 
into an appropriate context, they are compared in Q}able 4 
with corresponding results from tliree earlier experiments 
with the GISS model, as reported by Spar et al. (1976) 
and Spar (1977). Table 4 shov/s the rms errors; and SI 
shill scores, for both the model (F) and climatology (M), 
of sea-level pressures and 500 mb heights, as v/ell as 
rms errors of 850 mb temperatures, over the Northern 
Hemisphere for four months; January 1975, 1974, and 
1975, and February 1976. (The February F values are 





O’able 4. 


Soa-leve^ l>r ensure 
500 mb (m) 

I ;; ' 

1 850 mb (L’empQ 3 ?a 1 ;ure 

'Sea-level Pressure 
S 300 mb Height 



Summary of rrao errors and SI skili scores for 
foiu? GlSS-modoi simulations of monthly mean 
sea-level pressure, 500 mb heigh4, and 850 mb 
temporatare over the Northern Heiiii sphere* 

P denotes the model’ simulation and M repre- 
oents a ''forecast'' of climatology. 


Jan. 

1973 

Jan, 1974 

’4an. 

1975 

Peb. 

1976 

F 

M 

■ P M 

P 

M 

P 

M 

A. Bms 

i errors 





(mb) 10,0 

8.7 

8.6 9.2 

'5.3 

6.6 

8.8 

6.5 

72 

94 

80 108 

62 

82 

81 

99 

(K) 4.1 

4.3 

"4.7 5.1 

■' /:■' X r: 

4.1 

4.5 ; 

- 4.4 

4.6 

B. ^ 

scores 

/ .. 

-■ 





81 

81 

79. 89 

64 

; 73 

75; 

81 

■.} 

45 

55 

55 fSO "4 

'42^. 

52 

43 

53 


averages of PI and P2, ) In^terms of sea-level pressure, 
the model exliibits little 6r no skill relative to clima- 
tology, In two of the four months (including Pehruary 
1976) the model rms errors in sea-level' pressure exceed 
those for climatology. \Vhile the model •SI scores for . 
sea-level pressure are all equal to or less than those 
for climatology, their magnitudes are generally so large 
(except possibly in 1975) as to indicate no skill in the 
reproduction of the pressure patterns. On the other 
hand, the model does exliibit small but consistent skill 
relative to clijnatology in its. simulation of the montlily 
mean 500 mb height and 850 mb temperature fields. 

The global forecast (top) and observed (bottom) 
monthly mean sea-level pressure fields for Pebruary 1976 
are shown in Pigure 5, Major defects in the simulation 
are found in the Kortliem Hemisphere sub-tropinal highs 
(too wealc), the Icelandic Ipv/ (too far south), ajid the 

r.)i 

Eurasian high (too weak), ,, /as well as in the Southern 
Hemisphere- where the sub-polar lov^ pressure belt is 

/f 

underestimated. In Pigure 4, showing the forecast (top) 
and observed (bojbtom) 500, mb height fields for Pebi-uary 
1976, the simulation errors appear most serious over 

li-. ■' ' ■ " ’ 

Europe and the eastern North Atlantic Ocean, where the 
model fails to reproduce a large amplitude ridge. The 
failure of the model to simulate adequately the anomalously 
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war*a regime over North America in Pehmary 1976 is illus- 
trated in Figure 5# which shows the error (fo“i??ecast r 


observed) in the 850 mb temperature pattern over that 
region. Here it can’ be seen that negative errors as 

k 

large as - 10 K are found near the Great Lalces, and errors 
as large as t 12 IC appear over Alaska, Apparently the 
model is not yet capable of reproducing such a major 
climatic anomaly from the given initial conditions. 

(Ihe large positive error over Baff inland is yet another’ 
example of a defective simulation of the temperature 
field. )_ _ 


4 . S ea^surface I'eiaperature Ajiomaly IDxnerime nt ( P3 ) 

Interactions betv/een atmosphere a:ad ocean, and tiie 
influence of anomalous sea-surface temperatures (SST^s) 
on the behavior Of the atmosx)here, especially over long 
periods of time, have been the subjects of innumerable 

r-'.. 

studies,, speculations, midj iblications . As' part of 
the present investigation, a third, monthly meali forecast 
(P5) was computed for Pebriiary 1976, v/ith observed SST’s 
used in place of the climatological mean February values. 
Tlie purpose of this experiment was to measure the in- 
fluence of the SST anomaly field on the predicted monthly 
mean state of the atmosphere, and to determine if a fore- 
loiowlodge of the global SSI field would result in an 

' ^ ^ fr 
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improved oimulation, 

lUhe effect of a variable ocean surface on the evolu- 
tion of the atmosphere should bo studied via an inter- 
active oceEin-atmosphere model j but, as such a model is 
not yet available in practical form, an alternative pro- 
cedure is to assume that the 'SS'J? field is predictable, 
ahd'Lto use the observed sea temperatures as a surface 
boundary condition, ©le forecast designated v/as com- 
puted on this principle. 

In an earlier' experiment (Spar et al., 1976), it 
v/as found that insertion of observed daily SSI’s during 
a forecast run had no detectable beneficial effect on a 
monthly mean simulation. Ihe present SST anomaly ex- 
periment, ]?2., differs from the earlier one in that the 
daily observed SSI’s v/ere averaged over the month, and 
the monthly mean observed SSI field for Pebruary 1976 was 
then used in place of the climatological field as a fixed 
surface boundary condition, (Hie use of monthly mean ,, 
SSI’s instead of daily updated values v/as expected;'to „i., 
reduce the* ’’noise" effects of the data.) 

Ihe daily SSI data were provided. by the National 
Environmental Satellite Service (NESS ), and are based 
largely on satellite scanning radiometer measurements 
v/ith quality control through use of ship observations- ' ' 
(Brower et al., 1976), Ihese daily fine-mesli"NESS data 
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wore interpolatod to tlio 4° x 5° G-ISS grid, and avex-agod 
over the month, IQio deviation of the monthly mean Poh- 
ruary 1976 SSffi field from t)io climatologioal pattern for 
the month io illuistratod in Pigure 6, wliiclT-pliowa the SSU? 
anomaly only over tiio northv/eat quaGremt of the eartli, 

!Oic SSf hnomaJ.y iaotherms, dra\m for an interval of 1° 0, 
reveal oevoral poolo of relatively warm water, with 
raarima of + 3° 0 in the eaotorn faeific' and + 4° 0 in the 
weatern Atlailtic. IDhe global map (not shown) e^iibits 
smaller anomalies (up to -i- 0) in the western Uorth 

Pacific, and a bmid of positive anomalies (up to 4 5^ G) 
in the subtropical latitudes of the Southern Hemispliere, 
!DJie impact of SSf anomalies on the atmosphorie 
mo')tlily mean simulation is indicated numeric allj'’ in 
Sable 9 for the sariie regions as in SPablos 2 ‘and' in 
terms of rms differenoes and SI comparison scores be-- 
tween fox'ecasts P3 and P2, for sea-level 'pressure , 500 mb 
height, and 850 mb temperature. Comparing 2?ibles 5 and 

- i' 

2, it is apparent that the rms impact of the SSI' anomaly 
on the oee>lovel pX’essure field is considerably larger' 

j; 

than the rms replication error over all regions.,. How- 
ever, tile impact on the sea-level pressure pattern as 
represented by 'the Si scores is only slightly greater 
than tlic erroi' ^f 'Replication, and is therefore of ques- 
tionable significance. !J3ie 500 mb data in Table 5 
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Fig. 6. Monthly nean sea-surface temperature anomaly (observed-cl.’natology) for 
February 1976 over the northvest quadrant. 1®C isotherms. 
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Globe 
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Nortliern Hemisphere 


49 

39 

26 

1.6 
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Tropics 

5.0 

40 
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' E, (Pacific-U. Si 

6.2 


54 
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North America ; ' 
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56 

52 
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United States 

6,0 

54 
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Europe 

-4.7 

51 

37 

32 
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similarly exhibit a relatively largo rmo impact (compared 
with the replication results in I’ahle 2), but no signlfi- 
cant effect on the contour patterif ;as represented by the 

' n \ > 

SI scores, whicli cure about equally small in Tables 5 and 

-^T) jl \ 

2. Modest impact is also indicated bn the 850 mb tom- 

peratures by the rins differences in Table 5, which are 
only slightly greater than the reirlication errors in ' ' '' 

Table 2^. 

Miether or not the small impact of the SST anoma- 
lies on the atmospheric simulation v/as beneficial may be 
judged from the rms errors and Sl^ skill scores for fore- 
cast PJ, shown in Table S, which should be compared v/ith 
the corresponding error statistics for H and 1*2 in , 
Table'?. In general, a comparison of Tables 6 and ? re- 
veals that the use of '‘observed” SST*s in place of 
climatological values did reduce the simulation errors 
over' the globe and Northern Hemisphere, but not neces- 
sarily over the smaller sub-regions. IVhen evaluated 
over the globe and Northern Hemisphere, all the errors of 
shown in Table 6 are smaller than the corresponding 
errors for N1 and "S2 aho\m in~ Table 5, Furthermore, 
the rms errors of 500 mb height ax'e /ionsistently smaller 
for ,F3 than for FI and F2 over all regions (although the 
same is not true of the SI seCi'es). However, when one 
examines the sub-regions, 3 through 7, it appears that • 
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Table 6 « Root-mean- square (rms) errors ajid SI skill 
scores for UTebruary 1976 forecast 
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Sea-lerel 

Rressuro 


500 mb 
Heifflit 


850 mb i 
Temperature 


Region 

rms (mb) 

SI 

rms (m) 


' fms 

1. Globe 

5.8 

60 

51 

"38 


2, Northern Heitiisphere 

7.2 

72 

66 

fi 

41 

4.3 

3. Tropics 

3.4 

57 

12 

59 


4. E. Pacific-U. S. 

5.7 




5.8 

5» Nor th America - 

9.5 

86 ' 

86. 

52 



S. United States 

6.1 

' 1^, 

96 

84 

28'^ 

6.6 

7. Europe 

14.1 

108 

177 

68 
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the T5 simulation errors are actually worse than those 
for I*! ana. P2 in many instances, notably over the eastern 
Pacific and North Ainerica. (Chis is particularly true for 
the sea-level pressures andi 8!50 mb temperatures. 

Some of the synoptic effects of the SSI anomaly 
field are illustrated in■I'ig^lreB 7 and 8, which display 
the global simulations for of the monthly mean sea- 
level pressure and 500 mb height fields respectively. 
Compared vdth Pigure Figure 7 shov/s stronger and more 
realistic subtropical highs, but a v/e alter _a‘r'id less real- 
istic North Pacific low. 'Bie 500 mb simulation for- P5, 
shov/n in Figure 8, appears hardly distinguishable from 
the forecast in Figure 4, indicating no obvidiis major 
impact of the SSf anomaly, fhe effect of the SST anomaly 
on the 650 mb temperature field is illustrated in Figure 
9, which shov/s the difference between the monthly mean 
850 mb temperatures generated by F3 and by F2 (F5 - F^2) 
over the northv/est quadrant. The positive differences 
over Alaska and northwest Canada, v;hen compared vdth the 
error field in Figure 5» do indicate an improvement in 
the temperature simulation over those areas. ^However, 
thdso local effects are not reflected in tKe'Wgional 
error statistics of Table 6. 
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Fig. 7. Forecast aonthly aean sea-level pressure field for February 1976 co-outed 
vlth observed SST's (?3), 4 =b isobars. 


^nvnf) tinoff 




Fig. 8. Forecast r.onthly r.ean 500 r.b height field for February 1976 computed with 
obser%'ed SST's (F3). 100 M contours. 




Difference (F3-F2) between r.ean February 1976 850 cb temperatures computed 
with observed (F3) and climatological (F2) monthly mean SST's. I'K isotherms. 



5. C6nc3.usions 
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3!he results of the replication experiment indicate 

I'l -i 

w ' . , 

that small inherent uncertainties may exist in a monthly 

ra'ean simulation of the atmosphere computed v/ith. a global 
GCM. Qliese may be duo to minor computational difficul- 
ties > such as. interruptions and restarts , v/hich may 
prevent the exact duplication of the chain of computa- 

■■ S') 

tions leading to the final output. In the GIS3 model 
experiment, the rms errors of replication over the 
Northern Hemisphere were found to be about. 2 mb, 20 m, 
and 1 K for sea-level pressure, 500 mb height, arid 850 
mb temperatu^l respectively. In terms of SI compari- 
son scores, the replication errors of sea-level pressure 
and'^ 500 mb height over the Northern Hemisphere v/ere found 
to be about 40 and 20, respectively.. These results may 
be viev^ed as provisional estimates of the computational 
*'noise" level of the model, below which experimental 
••signals" cannot be reliably detected. 

0]he replication errors are small compared with the 
errors generated by the model in its effort to simulate 
the real monthly mean atmosphere. Measured in terms of 
rms errors and SI shill scores, and compared with the 
standard of climatology, the model simulation exhibits 
" no skill in reproducing the monthly mean sea-level 
pressure field. However, the model does show some small 



but coneistent skill relative to climatology in its simu- 
lation of the global and hemispheric fields of 500 mb 
height and 850 mb temperature. Unfortunately, this 
skill does not extend to the correct prediction of regional 
climatic anomalies such as the extreme v/armth over North 
America during Pebruary 197^. 

The use of observed montlily mean sea-surface tempera- 
tures (SST’s) as lower bomidary conditions in place of 
climatological values yielded ambiguous results. Tne 
impact of tiie SST anomalies on the monthly mean model 
-aimulations v/as relatively small, in general, and not 
consistently _beneficial. \^ile an improvement in the 
simulation scores over the globe and Northern Hemisphere 
did resuit from the use of observed SST* I, i this benefi- 
cial imijact v/as not consistently reflected eithex' in^-“. . 

il '' ' i'i 

the regional scores or jji ‘bhe monthly mean synoptlb mapl 5 . 

Purther improvements in the G-ISS model, as v;ell as 

'I 

. in the specification of surface boundary conditions 
(e.g., snow cover, soil moisture), are obviously desir- 
able. TOiether this vdll result in better monthly mean 
^ simulations remains to be seen, 

/f' 
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PIG-URES 


1. Hean meridional profiloo of vertically averaged mean 
zonal v;ind£3 for February 197o. PI and P2 repreaent 
the model simulations, and 0 is the observed profile. 

2. Vertical profiles of mean zonal Winds at latitude 
50°H for Pobruary 1976, PI and P2 denote model 
simulations. 0 is the observed profile, 

5. Porecast (top) suid observed (bottom) mean sea~level 
pressure fields for Pebruary 1976, 4 mb isobars. 

4. Porecast (top) and observed (bottom) mean 500 mb 

height fields for Pebruary 1976, 100 m contours. 

5. Error (forecast - observed) of 850 mb temperatia’e 
simulation over the northwest qtuadrant, 2 iC 
isotherms , 

6,, Monthly moan sea-surface temperature anomaly 
(observed - climatology) for Pobruary 1976 over 
the northv/ost q^uadrant, 1^ G isotherms. 

7. Porecast monthly mean sea-level pressure field for 
Pebruary 1976 computed with obsei'ved SST’s (P5)* 

4 mb isobars. 

8. Porecast monthly mean 500 mb height field for 
Pebruary 1976 computed v/ith observed SSI’s (P5) * 

100 ra contours. 

9. difference (P5 - P2) between mean Pebruary 1976 
850 mb temperatures computed v/ith observed (P5) 
and climatological (P2) montlily mean SST’s. 

1 K isotherms. 
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